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TECHNICAL MEMORANDUM NO. 542. 

CONTRIBUTION TO THE AILERON THEORY.* 
By A. Betz and E. Petersen. 

In an attempt to treat theoretically the effect of ailer- 
ons, difficulty arises because an aileron may begin at any 
point of the wing. Since the deflection of an aileron has the 
same effect on the wing as increasing or decreasing the angle 
of attack, a wing with aileron in action behaves like a wing 
with irregularly varying angle of attack. From the wing theory 
it is known, however , that the lift at such a point with irreg- 
ularly varying angle of attack does not vary irregularly. Hence 
the question arises as to how the transition af the lift dis- 
tribution proceeds at such a point, since the effect of the 
aileron (i.e., the moment generated about the longitudinal axis) 
depends largely on this distribution. 

In order to answer this question regarding the lift distri- 
bution during irregular variations in the angle of attack at 
first independently of other influences, especially those of the 
wing tips, we have taken as the basis of the following theoret- 
ical discussion a wing of infinite span and constant chord which 
exhibits at one point an irregular variation in the angle of 

*"Zur Theorie der Quer ruder" from Zeitschrift fur angewandte 
Mathematik und Mechanik, Volume VIII, 1938, pp. 353-257. 
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attack.* As regards the mathematical treatment, we will first 
consider a wing with periodically recurring irregular angle of 
attack (upper part of Fig. l). Ultimately we can let the period 
extend to infinity and then obtain the desired result for an 
infinitely long wing with a single point of irregular variation 
in the angle of attack. The treatment of a periodically varia- 
ble wing offers the advantage that the functions involved can 
be expressed in a Fourier series, which gives especially simple 
relations in the. present case. 

In order to express the lift distribution, we will seek the 
circulation V in terms of the distance x from the point of 
disturbance. Between the circulation V and the lift per unit 

length there is known to be the relation 

dx' 

M=pvr (i) 

dx 

in which p is the air density and v the flight speed. Ac- 
cordingly the lift coefficient at the given point, is 

*The application of the results to wings of finite span is disrr 
cussed by E. Petersohn, "Theoretische und experimentelle Unter- 
suchungen der unter Einwirkung von Querrudern an Tragflugeln auf- 
tretenden Moment e, 11 Luftf ahrtf orschung,' Vol. II, No. 2. 

Another treatment, based on an elliptical wing, was accorded 
the aileron problem by Dr. Max M. Munk (N.A.C.A. Technical Report 
No.. 191: "Elements of the Wing Section Theory and of the Wing 
Theory," 1924). 

While the present article was in press, another article, "The- 
oretische Unter suchungen uber die Querruderwirkung beim Trag- 
flugel," by C. Wieselsberger , appeared on this subject (Report 
No. 30 of the Aeronautical Research Institute, Tokyo Imperial 
University). In this article the lift distribution over a wing 
is approximately represented by a finite series of only eight 
terms. 
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c„ = 



d a . sr (2) 



a f v* t dx vt 

in which t represents the wing chord. 

The lift coefficient of a wing section or profile in an un- 
disturbed two-dimensional flow, can, with sufficient accuracy, be 
assumed to be # a linear function of the angle of attack a .* 

c a = c (a - a 0 ) (3) . 

d Co 



Thereby c = 



da 



a characteristic constant of the wing section. For flat plates 

the theoretical value is c = 2tt ; for thicker wing sections it 

is -somewhat greater* The m actual values are somewhat smaller 

than the theoretical. 

From equations (2) and (3) we obtain, the relation between 

T and a . 

r = c ^ (a - a 0 ) (5) 

whete a 0 is the angle of attack at which c a = 0*- The angle 
of attack of the wing may vary irregularly from a x to ct 2 (upper 
part of Fig. l). The circulations corresponding to these angles 
of attack in undisturbed flow (i.e., for an infinitely long wing 
with constant angle of attack) are then 

r\ = c ^ (a, - a 0 ) (6) 

and 

T 2 = c ^ (a a - oo) ' (7). 

♦Naturally this does not hold - true in the vicinity of the bur- 
ble point or after the flow has separated from the wing. 
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For reasons of symmetry a mean circulation ^ will 

prevail at the point of irregularity. The circulation from 
there on will approach asymptotically the value V x on one side 
and T 2 on the other side. We can therefore write 

r = - 3 ? + 1 3 a e ( 8 ) 

in which c is a temporarily unknown function of x. Our task 
is to determine the function e (x). 

The process of calculation is as follows. We develop a 
in a Fourier series and put T likewise in the* form of a Fourier 
series with temporarily unknown coefficients* From this distri- 
bution of r we can calculate, by the well-known wing theory 
method, the vertical induced velocities w on the wing, which 
alter the effective angle of attack by the amount 

A a = - | . (9) 

so that the effective angle of attack is 

, w 
a' = a - - . 

The circulation at every point x of the wing is calculated 
from this effective angle of attack according to equation (5). 
Since all functions are represented in the form of Fourier series 
the circulation distribution thus calculated is in the form of 
a Fourier series. The still undetermined coefficients of this 
series can be found by comparing the calculated circulation dis- 
tribution with that originally assumed. 
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The series for an irregularly varying angle of attack is 

a = + izfs i(sin ^ + J sin 3 + J sin 5 *£E + . . . ) 

o 2tt to 1.5 t 

do) 

(Of. Hutte, 25th edition, Volume I, p. 169.) For the dis- 
tribution of P we write 

P +r p R -r 2 , 2nx 2nx , c 2w . 

r = 1 0 g + 1 2 sm — ^ — + 03 sm3 — — + a 5 S m5 — — + ...) 
2 2 I I I 

From the distribution of V and according to the well-known cal- 
culation method of the wing theory the induced velocity w be- 
comes 

w= 1 /"ai _2= — dx (12) 

at a point on the wing x x distant from the point of disturb- 
ance. The summation of T according to equation (ll)* gives 



w = 



Fl g ^ *(2n + 1) £(2n + 1) sin (2n + 1)^5 (13) 



Since this calculation naturally holds good for any distance x x 
and not simply for a certain fixed distance, the subscript 1 
may be omitted and equation (13) would then represent in general 
the relation between the induced velocity w and the distance 
x from the point of disturbance. 

We may express the effective angle of attack 04 = a - w/v 
as a function of x and from it calculate the circulation V 



*L. Prandtl, "Tragflugeltheorie 1 ,; Fart I, Vier Abhandlungen zur 
Kydrodynamik und Aerodynamic, Gottingen, 1927, published by J. 
Springer, Berlin. _Under No. 14 it is shown that a circulation 
distribution T = T cos u- x gives an induced velocity 



w = J±fcos|i x = T . 
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P _ n vt / w\ _ vt 



a i + <h 



a n + 



a - a„ 4 

^1 2 



_ £ — sin (2n + 1) ^ 
tt q 2n + 1 t 



r. _ n, 



- c I -±-^-2- g cc (2n+l) IL (2n + 1) sin (2n + 1) ^J* (14) 

If we consider that, according to equations (6) and (7), 

vt / a i + a 2 _ \ _. c i + r 2 _ vt «i - °s _ r i - r 2 



c ^ f 
C 2 \ 



- a 0 ' = 



and c — 
3 



we obtain, "by comparing the last equation with the original 
equation for V (equation 11), the following relation 



r - r 

X l 1 2 



I a (2n+l) sin < 2n + D *t 



2ttx 



r i - r 2 00 



^(^Vry - t a (3n + D sln < 2n+1 > ¥ (15 > 



Since the coefficients of the corresponding terms of the two 
Fourier series must be the same, we obtain, for the coefficients 
a(2n+l) of the original summations for V , the expression 

4 ctn 



= 4 1 

a (2n+l) tt (2n + 1) 1 + ctn (2n + ]_) 

4 I 



(16) 



For the desired function e in equation (8) we therefore obtain 



0 tt (3 n + i) l + ctn (2n + x ) 



sin (2n + 1) 



2nx 



(1.7). 



■ -j i 



i 
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Thus we have solved the problem for periodical alternations 
in the angle of attack. In order to adapt the results to the 
case of a single point of disturbance, we must let the period 
I extend to infinity. For very large values of I and small 
values of x all the terms in the above series having small 
values of n approach zero as a limit. For large values of n, 
however, since n eund n + 1 differ but little, we can replace 
Z by an integral by introducing a uniformly varying quantity 
\ in place of the whole numbers n, so that 2n + 1 = 2 \ . 
The series (equation 17) then becomes 

4 T sin 2 \ p, " . 

tt 5 2 \ (1 + 2 v X) . UO; 



where, for brevity, we put 



I 

and 

v = 



\X = (19) 



C t TT (20) 



4 I 

This integral can be reduced to the well-known functions* sine 
integral 

S i l = /t sin z d z (2i) 
6 . z 

and cosine integral 

0 i | = £ d z (22) 

By partial fractional resolution the integral of equation 
(18) can be transformed into 



*E. Jahnke and F. Emde, "Funktionentaf eln rait Formeln undKurven," 
Leipzig, B. G. Teubner, 1923. 



i 
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r sin 2 X u dX = f sinUji dX _ f g sin 2- X jx ^ 
0 2X(l + 2vX) o 2X e 1 + 2 u X 

By the introduction of 2 X (i = z the first integral of the 
right-hand member of the above equation "becomes 



x / sin_z dz = ± n 
<d c z 4 

and by the introduction of ^ (l + 2 v \) = y the second inte- 
gral becomes 

I f /v (cos ^ sin y - sin B 00s y ) f . | cos £ (± J - Si H) + 

+ i sin 14 Gi H .* 
2 u u 

By the further introduction of *j = ff (equations 19 and 20) 
we obtain 



e = g 

TT 



± 3 ; cos ^ + cos ^ Si ^ - sin if Oi ff 
2 2 ct ct ct ct ct m 



= (± 1 - 2 sin 8x ci 8x^ „ cos 8x /+ i _ g si §x\ 
\ n ct ct/ ct V tt ct/ 

The behavior of the function e for positive values of x is 
shown in Figure 2. Negative values of x give the same curve 
but with the opposite sign. For large, values of x the function 
e(x) can be represented by the semiconvergent series 

, = ± ! _ J_ (i - 

• TT Z \ 



_2_ fl - 21 + 41 ^ 

2 Z / 



into which z = ^ has been introduced for. brevity. For small 

ct 

values of x the function is represented by the expression , 

*The positive sign corresponds to positive \i and x; the nega- 
tive sign to negative |i and x. 
**The series can be used only so long as the terms decrease. 



{■•• . • 

.... . ■;*■ 
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e =| (1 -C - In z) z 
(for z<^ 1), in which G is the Suler constant = 0,577. 



Translation by- 
National Advisory Committee 
for Aeronautics. 
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